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Abstract 

We compare two non-perturbative techniques for calculating the single- 
particle Green's function of interacting Fermi systems with dominant forward 
scattering: our recently developed functional integral approach to bosoniza- 
tion in arbitrary dimensions, and the eikonal expansion. In both methods the 
Green's function is first calculated for a fixed configuration of a background 
field, and then averaged with respect to a suitably defined effective action. We 
show that, after linearization of the energy dispersion at the Fermi surface, 
both methods yield for Fermi liquids exactly the same non-perturbative ex- 
pression for the quasi-particle residue. However, in the case of non-Fermi liq- 
uid behavior the low-energy behavior of the Green's function predicted by the 
eikonal method can be erroneous. In particular, for the Tomonaga-Luttinger 
model the eikonal method neither reproduces the correct scaling behavior of 
the spectral function, nor predicts the correct location of its singularities. 

PACS numbers: 72.10.Di, 05.30.Fk, 11.15-q,67.20+k 
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Typeset using REVTeK 



I. INTRODUCTION 



Recently a number of authors have developed generalizations of the bosonization method 
to arbitrary dimensions |[l|-[TO[l. This approach is most suitable in interacting Fermi systems 
where the Fourier transform of the interaction is dominated by momentum transfers small 
compared with the Fermi momentum kp- In this case the single-particle Green's function 
satisfies a non-trivial asymptotic Ward-identity [|ll|,^ , which opens the way for a controlled 
summation of the entire perturbation series. Bosonization is one possible (and perhaps the 
most efficient) way to explicitly carry out this summation. Interactions that exhibit in 
Fourier space a strong maximum or even a singularity for small momentum transfers are 
long-range in real space. The experimentally perhaps most relevant interaction of this type is 
the effective current-current interaction due to the coupling between fermions and transverse 
gauge fields. Such a coupling appears in models for strongly correlated electron systems 
TB| . Moreover, effective current-current interactions are also generated by fiuctuations of 



the statistical gauge field in the Chern-Simons theory for quantum Hall systems |14|. In 



Refs. [p|,p!5| the gauge field problem has been studied via higher-dimensional bosonization. 



However, this approach has been criticized |]T6|,|T2|, because higher- dimensional bosonization 



relies (at least in its simplest form) on the linearization of the energy dispersion at the Fermi 
surface, and there exists evidence that in the gauge field problem the curvature of the Fermi 
surface cannot be neglected |16|JT^ . 

An alternative non-perturbative approach, which does not rely on the linearization of 
the energy dispersion, is the eikonal method. This approach has been developed in the 
60 's to obtain non-perturbative results for the single-particle Green's function of field the- 
ories where the fermionic degrees of freedom are linearly coupled to another bosonic quan- 
tum field |[T7|-pO|. More recently, Khveshchenko and Stamp |^ have used this method 
to study the above mentioned gauge field problem in two dimensions. The strategy is to 
calculate first the Green's function Q{r,r' ,t,t'; {A}) of the fermions for fixed configura- 
tion of the gauge field A, and then average the result with respect to the effective gauge 



field action S'e//{A} to obtain the translationally invariant fermionic Green's function [22 



G{r — r',t — t') = {Q{r,T' ,t,t'; {A}))^ To leading order in the eikonal expansion the 
averaging involves a trivial Gaussian integration, and gives rise to the usual Debye- Waller 
factor |]2T|. The functional integral approach to bosonization |2^J^-|TD|,|T3[ follows precisely 



the same strategy: The fermionic two-body interaction is first decoupled via a Hubbard- 
Stratonovich transformation involving a space- and time-dependent auxiliary field; then the 
fermionic Green's function is calculated for a fixed configuration of the auxiliary field, and 
finally the result is averaged with respect to an effective action S^ff- Although the func- 
tional bosonization approach pB|,|7|- pmp^ is formulated in Euclidean (i.e. imaginary) time 
while the eikonal expansion is by construction a Minkowski (i.e. real) time method, the 
similarity in both procedures suggests that, after proper analytic continuation, they should 
yield equivalent results for models with linearized energy dispersion. In the present paper 
we shall carefully examine this point, and show that in general this is not the case. 
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II. FUNCTIONAL BOSONIZATION 



In order to compare bosonization with the eikonal expansion, we shall in this section 
briefly describe our functional bosonization approach, focusing on the calculation of the 
single-particle Green's function. For a more detailed description, see Refs. |f^^TO|. To cal- 



culate Matsubara Green's function of the many-body system, we subdivide the degrees of 
freedom close to the Fermi surface into a finite number of boxes K"" labelled by an index 
a, and decouple the interaction with the help of a Hubbard- Stratonovich transformation 
involving auxiliary fields associated with the i^". For our purpose it is sufficient to consider 
conventional density- density interactions, which can be decoupled via scalar fields 0". The 
generalization to gauge fields is straightforward and can be found in Refs. |T0|JT5[] . After the 
standard transformations P,pi]| the Matsubara Green's function G{k) = (^(k, icDn) can be 



exactly written as ||2^ 



G{k) = jv{r}v{r}[G]kk^{[G]tk)^^^^ . (1) 

Here is an infinite matrix in momentum- and frequency space, with matrix elements 
given by the formal Dyson equation [G~'^]kk' = [Go^]fcfc' — [y]kk', where Gq is the non- 
interacting Matsubara Green's function matrix, [Gojfcfc' = ^kk'[i^n ~ + lA~^^ ^'^'^ 
generalized self-energy matrix V is [V]kk' = J2a^"i.^)^k-k'^ with = ^(p^. Here /3 is the 
inverse temperature, /i is the chemical potential, ek is the non-interacting energy dispersion, 
and the cutoff function 0"(k) is unity if k belongs to box i^", and vanishes otherwise. The 
normalized probability distribution ^{0°} is 

where the effective action for the 0"-field is of the form S'e//{0"} = S'2{0"} + Skin{<P°'}, with 

Q aa' 

Sk^n{r} = -Tr\n[l-GoV] . (4) 

Here is a matrix in the patch indices, with matrix elements [/^]""' = yfq"'j where f^"' are 
the usual Landau interaction parameters, and V is the volume of the system. The probability 
distribution V{(p°'} can be calculated perturbatively by expanding Eq.| in powers of the 0"- 
field. The validity of this expansion is controlled by the generalized closed loop theorem. 



which is discussed in detail in Refs. PJTO[|. At the level of the Gaussian approximation one 
obtains 

Sk^n{r} ^ ^EE^T' iQ)r-,r,' , (5) 

g aa' 

where the "patch" polarization is given by 
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nr'(9) = E [e"(k)e"'(k + ^)Go{k)G,{k + q) 

+e-'(k)e"(k-q)Go(A;)G'o(fc-g)] . (6) 

The leading correction to Eq.^ is given in Refs. P,[lO|- The calculation of the diagonal 
elements [G]kk in Eq-|I| is more difficult. Choosing the patches larger than the range of the 
interaction in momentum space, we may write [G]kk = Y.a^"'(X)[G°']kk, where (after 

shifting k = k" + q, k' = k° + q') the infinite matrix satisfies 



q' 

Here [G^iq]]'^ = icD„-r(q), with ^"(q) = ek.+q-/i. Defining l^°(r,r) = Eg e^^'i'^-'^-^^l/" 
and 

6;°(r,r',r,r') = -L ^ e^{q r-^„r)g-.(q'.r'-^„,.')[G'"]_,, ^ (g) 

qq' 

it is easy to see that Eq.|^ is equivalent with 

[-dr - r (Pr) - V^"(r, r)] g^ir, r', r, r') = 6{r - t')6*{t - r') , (9) 

where Pr = — iVr is the momentum operator, and 6*{t — t') = ^ En Eq.|^ 
together with the Kubo-Martin-Schwinger (KMS) boundary conditions 



6;"(r, r', r + /3, r') = ^^"(r, r', r, r' + /3) = -^^"(r, r', r, r') , (10) 
uniquely determine the Green's function. For linearized energy dispersion, ^"(q) ~ v° ■ q, 



Eq.|^ can be solved exactly by means of a trivial generalization of Schwinger's ansatz |^ 
The result is [0-0 

6;"(r, r', r, r') = G^(r - r', r - rOe*"^'-'^)-*"^'-''^') , (11) 

1 pi(q-r-iljnr) 

pi{c\-r-UJmT) 

$-(r,r)=Vr C . (13) 

Note that the bosonic Matsubara frequencies in Eq.p!3| insure that $"(r, t + /3) = $"(r, r), so 
that the ansatz |ll] manifestly satisfies the KMS boundary condition |10[ Gaussian averaging 
of Eq.|ll| with the effective action S'e//{0°} given in Eqs.^ and |^ yields 

(^"(r, r', r, r')),^^^, = G^(r - r', r - r')e'3"(-r',.-.') ^ (14) 

1 £RPA,a 

^) II -c°e(q--^-.^)l . (15) 
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Here /^f^"^'" is the screened interaction within random- phase approximation (RPA), 



RPA,a 



L i + no(g)f 



(16) 



where underhned quantities are matrices in the patch indices. For the special case of the 
one-dimensional Tomonaga-Luttinger model Eq.^ exactly reproduces the well-known 



bosonization result 23,10 



III. THE EIKONAL EXPANSION 



An obvious disadvantage of the procedure outlined in the previous section is that the 



Schwinger solution in Eqs.|TT]-|T3] is only valid for linearized energy dispersion. In general, the 
effective energy dispersion in box K"" is of the form ek^+q = + v"" ■ q + The eikonal 
method can in principle handle the quadratic (curvature) term in a systematic way. However, 
this method cannot be directly compared with our functional bosonization approach, because 
it is by construction a real time method. Let us briefly summarize the main steps. For 
simplicity we focus here on the retarded Green's function for fixed real time background field 
V^"(r,t), which is defined via the differential equation 



[r,t) g:Ar,T',t,t') 



6{r - r')6{t - t') 



(17) 



with the boundary condition that ^"et(r, r', t, t') = for t — t' < 0. In Sec.0 we shall discuss 
how to obtain for the leading term in the eikonal expansion the corresponding time- ordered 
Green's function, which is usually calculated in diagrammatic perturbation theory. The 
eikonal expansion is a systematic method for solving partial differential equations of the 
type |l^. This approach has been developed many years ago [p!7|-pO|. In this section we shall 
briefly outline how Eq.^ can be solved with this method, following mainly Ref. |l^ . 

It is instructive to study first the solution of Eq.|l^ without external potential, i.e. for 
V^"(r,t) = 0. Then the retarded propagator is given by [^j 



^ret,o( 



r't 



_ r dq r ,[q.(r-r')-o.{t-t')]/7" („ 



(18) 



with 



G'ret,o(q.^) 



1 



/o 



(19) 



As will be evident shortly, the representation of the propagator in terms of an integral over 
an auxiliary time t" plays an important role in the eikonal approach. Note, however, that 
we can exactly eliminate the auxiliary integral by substituting Eq.|19| into Eq.|l^ and doing 
first the frequency integration, using / l^e"*"^*^*"* ~* = 6{t — t' — t"). Then the t"- integral 
is trivial, and we obtain 



<^?et,o(r - r',t - t') 



-ie{t-t') J 



dci 



i[q-(r-r')-r(q){t-t') 



(20) 
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If we linearize the energy dispersion (i.e. set = oo, as in conventional bosonization) , 
then the q-integration can be performed exactly, with the result 



r\t-t') 



m 



oo 



(21) 



Let us now consider the case V^"(r, t) 7^ . The starting point is the transformation of 
Eq.|l^ to the mixed representation by Fourier transforming it with respect to the difference 
variables r — r' and t — t' [|T^. Defining 



J[ciir-r')-w{t-t')]^a 



(22) 



{2ttYJ 27r 
it is easy to see that Eq.|l^ reduces to 

- r (Pr) - V^"(r, t) + u;- r (q)] G^r, t; q, 0;) = 1 . (23) 
The formal solution of this equation can be written as |jl9[ 

QUr, t- q, .;) = -^ H dt"e^'"^-^''^^^^^'^^Y'^it"; r, t) , (24) 



where the auxiliary function y°(t"; r,t) satisfies the partial differential equation 

'tdt"+idt-C{Pr)-V''{r,t)]Y''{t";r,t)=0 , (25) 

with boundary condition y°(0;r,t) = 1. The correctness of Eq.^ is easily verified by 
applying the differential operator in Eq.^ to the right-hand side of Eq.^ and integrating by 
parts. Note also that without external potential the unique solution of Eq.^ with the correct 
boundary condition is V^it"; r, t) = 1, so that Eq.|4| reduces to the right-hand side of Eq.|T9|. 
It is convenient to parametrize the solution of Eq.|2^ in the form Y°'(t";v,t) = e^"^*"'""'*^. 
The equivalent differential equation for is then 

[td,>+td,-C{Pr)]F"{t";r,t) = V"{r,t)+ ^ ^'^^''^ , (26) 

with boundary condition F°'{0; r, t) = 0. This non-linear partial differential equation cannot 
be solved exactly. However, we can obtain the solution as expansion in powers of the external 



potential H. Substituting the ansatz F°(t";r,t) = E^=i ^°(^"; r, ^) into Eq.|| (where 
F° is by assumption of order (K")"), we have 



+ zdt - r(Pr)] F:{t"- r, t) = Kit"; r, t) 
with V{'{t";r,t) = V"{r,t) (independent oft"), and 



n 



1,2,.. 



n-l 



C(t";r,t) 



2m° 



5:[Pri^„-,(t";r,t)]-[PrF„%,(t";r,t)] 



n = 2,3,... 



(27) 



(28) 



n'=l 



Eq.|23 is easily solved by means of the Green's function of the differential operator on the 
left-hand side. The solution with the correct boundary condition can be written as 



F:{t"; r, t) = -t t" t;se*[*^'-«"(P^)l(*"-^)C(s; r, t) . 
Jo 



(29) 



This completes the formal solution of the problem. 
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IV. FORMAL SIMILARITIES BETWEEN BOSONIZATION AND THE 
EIKONAL EXPANSION FOR LINEARIZED ENERGY DISPERSION 



Suppose we truncate the eikonal expansion at the first order, F"(t";r, t) ^ F"(t";r, t). 
This approximation becomes exact in the hmit = oo (i.e. for hnearized energy disper- 
sion), because then all higher order potentials with n > 2 in Eq.p7| vanish identically. For 
finite m", we expect that the relevant small parameter which controls this approximation 
is proportional to the product of 1/m" times the typical strength of the effective screened 
interaction P,|lOl- The function F^(t"; r, t) is easily calculated by decomposing V^"(r, t) into 
its Fourier components. 



and using 



gi[iat-e(Pr)](i"-s)gi[q-r~<^t] 



gi[a;-e"{q)](i"-s)gi[q-r-<^t] 



After performing the s-integration in Eq.|2^ we obtain 



i^r(t";r,t) 



(27r) 



277 u- ^"(q) 



,i[<^-5-(q)]t" 



(30) 



(31) 



(32) 



The leading eikonal result for the retarded Green's function in a given external potential is 
then 



Jo 



(33) 
(34) 



We now show that for linearized energy dispersion (i.e. for m" = oo) Eqs.|3^ and ^ have 
exactly the same structure as the Schwinger solution given in Eqs.|TI|-|I3|. At the first sight 
this is not at all clear, because the eikonal result in Eqs.^ and ^ involves an additional 
integration over the auxiliary time t" . In the work of Khveshchenko and Stamp |^T| the t"- 
integration has been performed by means of the saddle-point method. However, to see the 
connection with bosonization, we eliminate the auxiliary variable t" from Eq.0. Therefore 
we simply repeat the manipulations leading to Eq.^ and first perform the w-integration in 
Eq.33. As before, this gives rise to a factor of 6{t — t' — t"). The t"-integration is then trivial, 
and we obtain 



G?etAr, r', t, t') = G^^Ur - r', t - t^e^^" 



(35) 



with Q(r — r',t — t') given in Eq.^. Still, the exponential factor in Eq.^ does not have 
the same structure as the corresponding factor in Eq.|TT[ However, for m° = oo the function 
G°gjQ(r — r',t — t') is proportional to 5*^'^)(r — r' — v"(t — t')) (see Eq.21), so that we may 
replace v"(i(: — t') ^ r — r' in the expression for F^{t — t'; r, t) in Eq.p5|, 
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,i[a;-r(q)](i-i')' 



J[(l-r-ujt] 



I _ giaj(i-i')-iq-v"(t-t')' 



1 _ gi[<^(t-t')-q-(r-r')] 

Then it is easy to see that Eq.^ can also be written as 

G?etAr, r', t, t') = G^,,o(r -r',t- t')e''?ir,t)-nir',t') 

with 



27r 



(27r 



- V" ■ q 



(36) 



(37) 



(38) 



The formal similarity to Eqs.|TT|-|T3| is now evident. Note, however, that the imaginary time 
finite-temperature Green's function discussed in Sec.|I| and the real time retarded Green's 
function discussed here are not related in a trivial way: only after Fourier transformation 
into the frequency domain there exists a simple relation via the spectral function. 



V. AVERAGING IN MINKOWSKI TIME: HOW GOOD IS THE EIKONAL 
EXPANSION FOR LINEARIZED ENERGY DISPERSION? 

In spite of the formal similarity between the real time eikonal result in Eq.^ and the 
imaginary time bosonization solution in Eq.^, after the averaging procedure the final ex- 
pressions for the Green's function of the many-body system are not equivalent. To average 
Eq.|3^, one should keep in mind that by construction functional averaging always generates 



time- ordered correlation functions |2^. Hence, in order to properly define the averaging 
procedure, we should replace the retarded Green's function in Eq.0 by the corresponding 
time-ordered Green's function before averaging with respect to the effective action Se/f of 



the background field. This is easily done by going back to Eq.|T9| and choosing the limit for 
the ^''-integration such that the integral for the time-ordered Green's function is convergent. 
With the notation = signer we have 

G-(q, Lu) = -—\ — = -z rft"e^t"[-r(q)+^o+..] _ (39) 

' - e(q) + Jo ^ ' 

Using the fact that functional averaging restores translational invariance in space and time, 
it is easy to see that after averaging the leading eikonal result for the time-ordered Green's 
function reads 

G?(q, u) = -I rft"e**"[-€'^(q)+*o+s.] /F^if'm) _ . (40) 

At the level of the Gaussian approximation the functional average in Eq.^ generates the 
usual Debye- Waller factor, with propagator given by the RPA-interaction Thus, 



Jo 



with 
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^^^^^-y (2^7 2^^ ^^'""^ [c.'-e(q')]^ ■ ^ ^ 

For simplicity we have assumed a patch- and frequency-independent bare interaction /q, so 
that the screened interaction can be identified with the usual RPA interaction, /'^^^(q, co) = 
/q[l + no(q, co')/q]~^, where no(q, a;) is the Lindhard function In Eq.^ we have called 



the integration variables q' and u' to avoid confusion with the external labels q and u in 



Eq.|41 



According to Eq.^ the spectral function — ^/mG"(q, u + iO"*") of our interacting many- 
body system can be obtained via a one- dimensional Fourier transformation. Furthermore, 
except for the trivial factor of s^, the spectral function depends exclusively on the combi- 
nation u — ^"(q). In particular, for linearized energy dispersion ^"(q) = v° ■ q (where the 
eikonal expansion truncates at the first order, so that Eq.^ is the exact solution of Eq.|T^ 
the spectral function is according to Eq.^ a combination of the variable a; — ■ q only. The 
well-known bosonization result for the Green's function of the Tomonaga-Luttinger model 
p7| (i.e. one-dimensional electrons with long-range interactions and exactly linear energy 
dispersion) shows that in one dimension this result cannot be correct |^! The spectral 
function of the Tomonaga-Luttinger model exhibits a more complicated dependence on q 
and uj One possibility to cure this problem might be to choose some effective u- or 



q-dependent cutoffs for the integration limits in Eq.52 However, such a procedure 

seems to require a certain amount of intuition and knowledge about the final result, and is 
certainly not satisfactory from the formal point of view. 

The obvious question is now whether the spectral function calculated from the real time 
eikonal expansion given in Eqs.^ and ^ agrees at least approximately with the spectral 



function calculated from the imaginary time bosonization result given in Eqs.l4 and 15. To 
clarify this point, it is useful to rewrite Eq.^ with the help of the dynamic structure factor 
which is related to the RPA interaction via [|1^] 

f^^^ic{,u') = U,-fi rd^S^^^ic^,uj)\^-^ + ^-] . (43) 

^ Jo lu — Uj' U + Uj'i 



Substituting this expression into Eq.^ it is obvious that we encounter poles on the real 
cu'-axis, which must be regularized by an appropriate deformation of the integration contour 
in the complex cj'-plane. The correct contour is easily determined from the requirement that 
the expansion of the factor e'^"^^ ^ in Eq.^ should reproduce the lowest order perturbation 



theory for the time-ordered Green's function |34]. It is easy to see that this implies that we 



use the dashed contour shown in Fig.|^, which describes time-ordering and corresponds to 



the regularization uj' ^ uj' + iO~''signco'' in Eq.E3|. The tu'-integral can then be carried out via 



contour integration. We obtain Qi{t") = R" — Sf{t"), with 

rfq ^2 r , S^^^{q,uj) 



f.a_ f dq , (q,^j 



(27r)'^''''yo (cu+lv^-ql) 
^t" f dq ,nPAr__^ f dq ,2 r^^. ^''''^(q,^) 



{271^ J {2tiY-"'Jo (cu+|v°-q|)2 

J e"^"^*" [(cj2 + (v" ■ q)2) cos(v" • qt") + 2icjv" ■ qsin(v" ■ qt")] - 2cj|v" ■ q| 1 
M (^-|v-q|)2 1 • ^''^ 
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It is easy to show that the term i?" agrees exactly with the zero-temperature hmit of the 



space- and time-independent contribution to Q"{r,T) in Eq.|T5 



1 £RPA,a 

^? = ii- JL- t^Et— ^^t:^ • (46) 



Because in a Fermi hquid the quantity e^i can be identified with the quasi-particle residue 
PIP^ , it is clear that for Fermi liquids the real time eikonal approach and functional bosoniza- 
tion yield exactly the same non-perturbative result for the quasi-particle residue. The proper 
deformation of the contour into the complex tu'-plane was crucial to obtain this result. If we 
had directly averaged the eikonal result for the retarded Green's function in Eq.^ (choosing 
the retarded contour shown in Fig.|l]), we would have obtained incorrectly i?" = 0. 

The time-dependent contribution S^{t") cannot be directly compared with a correspond- 
ing term in Eq.|T^. To examine possible discrepancies with bosonization, let us explicitly 
calculate Qi{t") for the spinless Tomonaga-Luttinger model with interaction parameters 
9i{^) = 5'2(ci) = /q- In this case the patch index a = ± labels the two Fermi points, and the 
dynamic structure factor is simply S^^^{q,uj) = Zq6{uj — tUq), with residue — — 



where 



and collective mode = ^JY^-^\^VF\q\. Here Fq = /q/(vrfi?) is the relevant dimensionless 
coupling. (It is understood that in d = 1 the vector q has only one component.) It is also 
easy to show ||10[ that in this case /'^^^(q, v"-q) = 0, so that the first term in Eq.^ vanishes. 
Because of the 5-function dynamic structure factor, the tu-integration in Eq.^ is trivial. To 
perform the remaining (one-dimensional) q-integral, we adopt the usual procedure |^ of 
replacing /q — ^ /o in the integrand and multiplying the integrand by an ultraviolet cutoff 
The q-integration is then elementary and we obtain after straightforward algebra 

QUt") = -(1 + |) In [1 + z{vF - VF)q,t"] - I In [1 + z{vf + VF)qct"] , (47) 

7 = ^ , ; (48) 

is the well-known anomalous dimension of the Tomonaga-Luttinger model, and "Vf = 
i/l + FqVf is the renormalized Fermi velocity. For u > the eikonal result for the spectral 
function of the spinless Tomonaga-Luttinger model can then be written as 

1 /.oo jt"[u;-v°-q+iO+] 

v4"(q,c^) = -i?e / dt" —J. ^ . (49) 

7o [l + ^(^y^-yp)q^t'f^^^^[l + t{vF + VF)qct"V^^ 

Note that the integrand vanishes as {t"y~^'^ for large t". Hence, for any finite interaction the 
integral in Eq.^ is convergent, and we may omit the convergence factor iO~^ in the exponent. 
Following Dzyaloshinksii and Larkin [^, we find that for small \uj — \"'- q| and < 7 < 1 to 
leading order A°'{ci,uj) ~ +y4"|u; — v"-q|'^, where and A" are finite complex numbers. 
This is in striking disagreement with the well-known bosonization result |30|, which predicts 
for < 7 < 1 a singularity of the type A°'{q,uj) oc |a; — vT+T^v" ■ q]'''""'^. Thus, the real 
time eikonal method neither reproduces the correct scaling behavior of the spectral function, 
nor does it predict the correct location of its singularities. 
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VI. CONCLUSION 



In this work we have compared the real time eikonal method with functional bosonization. 
In both methods the translationally invariant Green's function G of the many-body system 
is obtained by calculating first the Green's function ^ in a given background field, and then 
averaging the result with respect to a suitably defined Gaussian effective action, G = {Q). 
Although both methods produce very similar non-perturbative expressions for the final 
results for the functinal average {Q) are certainly not equivalent. In particular, in the case 
of the one- dimensional Tomonaga-Luttinger model the real time eikonal method does not 
even reproduce the correct scaling behavior of the Green's function. We suspect that this 
shortcoming of the eikonal approach is related to the averaging in Minkowski time, which 
cannot be interpreted as averaging over an ordinary probability distribution. On the other 
hand, the imaginary time averaging performed in functional bosonization is mathematically 
well defined, and leads to an expression for the Green's function which is exactly equivalent 
with the result obtained via operator bosonization 0,0. 

Because even at the level of the linearized theory the real time eikonal expansion can 
lead to incorrect results, we conclude that this method is not a good starting point for 
studying the effects associated with the curvature of the Fermi surface on the low-energy 
behavior of the Green's function of non-Fermi liquids. Inclusion of curvature effects is of 
particular importance in connection with the gauge field problem. It is therefore highly 
desirable to study the effect of the non-linear terms in the energy dispersion on the leading 
(linearized) bosonization expression for the Green's function entirely within the framework 
of the imaginary time approach described in Sec.|l|. Work along these lines is in progress. 
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FIGURES 
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FIG. 1. Possible integration contours in the complex w'-plane for the evaluation of Eq.42. The 
dashed line is the correct time-ordered contour. The solid line corresponds to choosing the retarded 
interaction in Eq.42. Both contours are infinitesimally close to the real axis. 
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